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Abstract 

This article offers a new approach for analysing the dynamic behaviour of distributions of 
charged particles in an electromagnetic field. After discussing the limitations inherent in the 
Lorentz-Dirac equation for a single point particle a simple model is proposed for a charged 
continuum interacting self-consistently with the Maxwell field in vacuo. The model is developed 
using intrinsic tensor field theory and exploits to the full the symmetry and light-cone struc- 
ture of Minkowski spacetime. This permits the construction of a regular stress-energy tensor 
whose vanishing divergence determines a system of non-linear partial differential equations for 
the velocity and self-fields of accelerated charge. Within this covariant framework a particular 
perturbation scheme is motivated by an exact class of solutions to this system describing the 
evolution of a charged fluid under the combined effects of both self and external electromag- 
netic fields. The scheme yields an asymptotic approximation in terms of inhomogeneous linear 
equations for the self-consistent Maxwell field, charge current and time-like velocity field of the 
charged fluid and is defined as an ultra-relativistic configuration. To facilitate comparisons with 
existing accounts of beam dynamics an appendix translates the tensor formulation of the pertur- 
bation scheme into the language involving electric and magnetic fields observed in a laboratory 
(inertial) frame. 

1 Introduction 

The intense international activity involved in probing the structure of matter on all scales, with 
particle beams and radiation, owes much to recent advances in accelerator science and technology. 
Developments in the production of high power laser radiation also offer new avenues for accelerator 
design and new diagnostic tools of relevance to medical science, engineering and the communications 
industry. A common theme in these developments is the interaction between charged particles and 
the electromagnetic field in domains where relativistic effects cannot be ignored. 

It is remarkable that many of the challenges that must be addressed in order to develop and control 
devices that accelerate charged particles have their origin in the interaction of particles with their own 
electromagnetic field. Despite the fact that the classical laws of electromagnetism were essentially 
formulated over a century and a half ago the subject of electromagnetic interactions with matter 
remains incomplete. This incompleteness has had concomitant effects on the development of quantum 
electrodynamics and renormalisation theory. At root, the difficulties reside in the recognition that 
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the quantum structure of matter at some scale is beyond observation. Furthermore, the classical 
description of the electron as a point particle leads to singularities in the Maxwell self-fields that 
inevitably create ambiguities in its interaction with the Maxwell field. The general consensus is 
that a useful domain of validity of the Lorentz-Dirac equation PQ 121 IS] , describing covariantly the 
radiation reaction on a point electron, can be accommodated by performing a "reduction of order" 
that effectively replaces the equation by a perturbative second order system for the particle world 
line. One must then decide whether higher order terms in this expansion should be maintained given 
the neglect of terms associated with the regular isat ion scheme. This approach is behind many of the 
successful applications of approximate radiation reaction dynamics, despite the somewhat delicate 
and unsatisfactory nature of the arguments that purport to support this approach. 

It appears that the analysis of systems involving radiating matter has progressed by a very 
successful symbiosis between experimental expediency and a variety of approximation schemes with 
proven effectiveness in different domains of validity. However, as schemes for accelerating charged 
particles in these devices become more complex and ambitious in their aims it is apparent that some 
existing theoretical models are inadequate for a proper understanding of new challenges. Such models 
are sometimes enshrined in large commercial computer codes that do not survive scientific scrutiny 
or have documentation that makes contact with established scientific literature difficult. As particle 
energies increase with the use of higher intensity laser fields this issue becomes critical and more 
reliable methods for accommodating radiation reaction must be found. 

A mathematically coherent formulation of a closed system of partial differential equations describ- 
ing the relativistic behaviour of charged matter with electromagnetic fields is inevitably non-linear 
and, in general, exact solutions satisfying causal boundary conditions are intractable. The fact that 
the effects of (coherent and incoherent) radiation on micrometre-sized charged bunches in a host 
of newly proposed advanced devices is a significant barrier to development indicates that existing 
approaches used to model such effects are inadequate and that a new look at the whole problem is 
timely. 

Many of the radiation problems alluded to above can be circumvented by working entirely in 
the language of fields rather than both particles and fields. This makes it possible to use intrin- 
sic tensor field theory and exploit to the full the power of differential geometry in a relativistic 
spacetime framework. The geometric field formulation offers a number of powerful computational 
advantages over existing formulations of beam dynamics. These include the maintenance of relativis- 
tic covariance (local Lorentz transformation between local (accelerating) frames are unnecessary), 
the use of curvilinear coordinates adapted to the geometry of particular problems (thereby facilitat- 
ing imposition of boundary conditions), exploitation of available symmetries and the formulation of 
coordinate-independent approximation schemes. 

In this article a new approach is offered for analysing the behaviour of distributions of charged 
particles in a coupled electromagnetic field environment. Guided by a simple model, an approxima- 
tion scheme designed for charged distributions containing ultra-relativistic particles 1 is explored. The 
exact equations of motion are derived from the vanishing divergence of a relativistic stress-energy 
tensor for matter and radiation. Before presenting the approximation scheme in detail we recall the 
genesis of the radiation reaction for a point charge, mention its shortcomings and argue that these 
may be overcome by treating a large collection of charged particles as a continuum. 

1 A novel definition of an "ultra-relativistic vector field" relying entirely on the light-cone structure of spacetime is 
given in Appendix IA.1I 
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1.1 Equations of motion for radiating charges 

In this paper the language of differential geometry is employed since it offers the most succinct 
mathematical framework for electrodynamics. Thus, all fields will be regarded as sections of tensor 
bundles over appropriate domains A4 of Minkowski spacetime endowed with a fixed metric g and a 
torsion- free metric compatible connection V. Sections of the tangent bundle over Ai will be denoted 
TTM. while sections of the bundle of exterior p-forms will be denoted ThPM.. For any vector field 
X denote by X the associated 1-form defined by X = g(X,—). The operator d will denote the 
exterior derivative and ix the contraction operator with respect to X. To facilitate comparisons 
with existing accounts of beam dynamics Appendix IA.2I translates the tensor formulation of the 
perturbation scheme into the 3-vector language involving electric and magnetic fields observed in a 
laboratory (inertial) frame. 

Due to its high symmetry, Minkowski spacetime admits a class of global charts that play a 
fundamental role in the following. A generic system of coordinates in one of these charts will be 
denoted {y^} G M 4 where \x ranges from to 3. In these coordinates the metric tensor takes the form 

g = -dy° (S>dy° + dy 1 ®dy l + dy 2 ®dy 2 + dy 3 ® dy 3 . (I) 

With y° > the field ^ is defined to be future-pointing. When convenient y° = ct, y 1 = x,y 2 = 

y,y 3 = z is written and the vacuum speed of light c is set to 1. The significance of this class of 

r - JL 

^ dyf 

£k»9 = (2) 



coordinates is that it offers a basis of symmetry generators {K^ = jfj;}: 



where Lx is the Lie derivative with respect to X. Such Killing vectors will be used to define energy- 
momentum densities and power-momentum fluxes associated with different field configurations on 
spacetime. 

If U is a domain of spacetime with boundary 

dU = Si + s 2 + n (3) 
for space-like hypersurfaces Ei, E 2 and J a closed regular 3- form on U (i.e. dj = 0) then 

J = [ dj=0. (4) 

dU JU 

Thus 

/ J= [ J- I J- (5) 

Such closed 3-forms imply conservation laws. 

In a vacuum the electric and magnetic fields are encoded into the 2-form T and the electric 
charge density and current are described by a source current 3-form j . The Maxwell field system on 
spacetime is 

df = (6) 

and 

d*g = -j (7) 

where * is the Hodge map associated with the spacetime metric tensor g and Q = e§T where e is 
the permittivity of the vacuum. 
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For any vector field W on spacetime and any Maxwell solution T define the "electromagnetic 
drive" 3-form 

r^ M) = Ui w T A*F-i w *TAT}. (8) 
If W is a Killing vector field then r^ EM ' ) is called a "Killing current" and 

dr^ M) = -i w TA 3 (9) 
e o 

follows from ( ffij) . and (jHJ. For each Killing vector field these equations establish a "local conser- 
vation equation" ( drw = 0) in a source free region (j = 0). 

If U is any time-like vector field (with g(U,U) = — c 2 ) one may write uniquely 

F = eAU + *(bAU) (10) 

where g(e, U) = g(b, U) = and 



I r ( EM ) = -e A g A U - l -{g{e, e) + g(b, 6)K(*1)- (H) 



Poynting 2 related the form e A b, in a source free region, to the local field energy transmitted 
normally across unit area per second (spatial field energy current or field power) and \{g{e, e)+g(b, b)} 
to the local electromagnetic field energy density. 

More precisely - L 4 EM) is the field energy associated with the space-like 3-chain £ and js J s i^r^ EM ' ) 
is the power flux across an oriented space-like 2-chain S. 

If X is a space-like Killing vector generating space-like translations along open integral curves 
then with the split: 

T^ M) =fi X AU + Q x , (12) 

where ijj\xx = and iuGx = 0, the Maxwell stress 2- form fix may be used to identify mechanical 
forces produced by a flow of field momentum with density 3-form Qx- In any local frame {X a e 
TTM.} with dual coframe {e b G TA 1 ^} the 16 functions T^ M ^ = ix b * r x^ M ^ where the frame 
indices a, b = 0, 1, 2, 3, may be used to construct the second-rank stress-energy tensor 

T (EM) = T (EM) e a ^ e b_ 

A standard application of these general notions is the derivation of the equation of motion of 
a charged fluid in an external electromagnetic field from the total stress-energy of the fluid and 
electromagnetic field. A thermodynamically inert (cold) fluid can be modelled with the stress-energy 
tensor 

T ^ = r ^MV®V (13) 
ce 

where M is a scalar number density field, m some constant with the dimensions of mass, V the unit 
time-like 4-velocity field of the fluid and g{V,V) = —1. Such a stress-energy tensor gives rise to a 
set of Killing currents 

4l = —9(V,K,)k(AfV) (14) 
ceo 



2 In a pre-relativistic context. 
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which are added to t^ M ^ to yield the total set of Killing currents for the interacting system. If one 

assumes that the electric current 3- form is j = q$N *V for some electric charge constant go an d that 
M is regular then the conservation laws 

dj = (15) 



yield the field equation of motion 



d(r^ + 4l) = (16) 

V V V = -^%vT. (17) 

m cr 

This equation must be solved consistently with the Maxwell equations ©, (J7|) to determine V, Af 
and T for prescribed initial and boundary conditions. 

If the source of T contains a 'point charge in an external field this approach must take into account 
that part of the electromagnetic field becomes singular on the world-line of the particle so the domain 
U above must exclude this line. Furthermore, one needs to postulate a stress-energy tensor for a stable 
point particle. Without further information application of Stokes' theorem then leaves an unknown 
contribution to the conservation law. However, experiment indicates that accelerating electrons 
(thought to be point particles) experience a reaction to their emitted radiation and that the power 
radiated is in good agreement with the covariant Larmor formula j3] . A number of methods have been 
devised in an attempt to accommodate this observed radiation reaction by adopting a regularisation 
procedure that permits the use of Stokes' theorem in a spacetime region that includes the particle 
world-line. Inevitably this requires some assumption of how the singular Coulombic stresses are 
compensated by stresses that are not of electromagnetic origin. Early approaches invoked delicate 
limiting processes and the use of both advanced as well as retarded solutions to Maxwell's equations 
Dirac [3] offered one of the simplest covariant regularisation schemes and the resulting equation 
of motion is known as the Lorentz-Dirac equation. The use of both advanced and retarded solutions 
is however unnecessary Hj . 

To highlight the assumptions in the single particle approach to radiation reaction an outline is 
given of how the Lorentz-Dirac equation can be derived using a coordinate system adapted to the 
light-cone structure of spacetime and the time-like world-line of a particle on a parametrised history 
given in the above chart as y^ = ^(u). For each value of u the unit tangent to the world- line is the 
4- velocity 

V = V»(u)^L (18) 

where V^{u) = ■ 

Assign to any event in spacetime not on the particle world-line, whose backward light-cone 
intersects this world-line at the event with parameter u, the set of coordinates {u, r, 6, <ft} with 
< r < oo, 0<#<7r,0<</><27T using the transformation 

y» = e(u)+ \ TT ^(M, (19) 

v{u,e,(p) 



where for any y° > : 
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5> fc (M)) 2 = i (21) 



k=l 



V = -g(V,L), g(V,V) = -l, g(L,L) = 0. (22) 

In these coordinates the particle world-line is the curve where 3 r = 0, (0, (f)) parametrise the unit 
2-sphere S 2 where {n k (9,(f))} are direction cosines and the metric tensor can be written: 



where 



g = - e u <g> e u + e 1 <g> e 1 + <g> e 2 + e A <g> e A (23) 



e° = Bdu- ^dr (24) 

e' = % (25) 
e 2 = ^-d# (26) 



e 3 = ^sinfld0 (27) 



^=\/ 1 + 2 ^. (28) 
Here, and in the following, denote d u f by / for any scalar field / and let T = Ca_T for any tensor 

du 

field T. 

Suppose the total system (electromagnetic field and charge) has a stress-energy tensor giving rise 
to a regular Killing current t k such that 

dr K = (29) 

in a domain U bounded by the light-like 3-chains E(u = w + Am), E(u = wo) and the time-like 
3-chain E(r = r ) for arbitrary positive constants Wo, A-u,r . Then 

/ rjf- / t k + f r K = 0. (30) 

In the limit Aw — > this may be written 

PP (u , r )du + PP (u ,r o )du = (31) 



where 



P£\u ,r ) = f [ i K {u o ,r,0A) (32) 
Jo Js 2 

pK\uo,r )= id u TK{u ,r ,6,<t)) (33) 
Js 2 



3 For V > with V time-like and future-pointing the set r > lies on the forward light-cone of the event 
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With K e {g^r} and = one has a similar balance of rates for each translational Killing 

vector in the basis. In terms of the "body" rate 

p(B) 



p(B) 
r a 



dy" 



and the "contact" rate 



there is a balance of 1-forms: 



P^ = P { ^dy^ 

dyV- 



P (B) K,r )+P (C) (uo,r ) =0. 



Next assume that the contact rate includes a part from the Killing currents tL \ i.e. 



p(C) _ p(C) , p 
— - r MECH ^ -H 



>(C) 
EM 



(34) 

(35) 
(36) 

(37) 



and furthermore that this part can be calculated from the retarded Lienard-Wiechert solution T to 
Maxwell's equations for an arbitrarily moving point charge qo in no external electromagnetic field. 
In these coordinates the solution is T = dA where the 1-form 



.4 



Qo V(u) 



47re r 

Integration over the 2-chain S 2 with r = r yields for (|33|) : 



EM 



4^ 



(3* 



(39) 



U=UQ 



where A = V^(u)-^: is the acceleration field. Clearly this 1-form is singular on the world-line where 
ro = 0. The first term however correctly accounts for the observed Larmor radiation rate of energy- 
momentum from an accelerating charge and is independent of tq. To proceed one must cancel the 
singular rate in Pg M from singular terms in the remaining rates in such a way that in the limit ro — > 
the resulting system of ordinary differential equations for arbitrary uq makes both mathematical and 
physical sense. One approach is to suppose that the remaining rates are determined in terms of 
scalar fields a(u, r), /3(u, r) and the vectors V, A such that 



+ 



a 



P {B \u,r ) + Pme CH (m, r ) = a(u, r )V{u, r ) + f3(u,r )A{u,r ). 
Inserting this in (|3^j) and applying the projection operator U v = 1 + V A iy one finds 

P = - 

With this value (}3l)j) yields 
A 



2_^_ 

3 4vre; 



A{A) 



(40) 



(41) 



a 




a 



A(A) 



+ 



4vreg 2r 



-V 



a + 



3 4vreg 



A(A) 



2jt 
3 4vre; 



+ 



a 



A(A) 



A 



(42) 



where A 



V^iu)-^. One way to simplify this and cancel the exposed electromagnetic singularity 



is to further assume that a(u, r ) 
one has: 



1 2 o 



1 



for some constant m so that in the limit r — > 



3 47re r 



(43) 
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or 

m c 2 A=-\-^i v (y A~A). (44) 
3 47re 

This is a manifestly covariant equation of motion even though the velocity and acceleration are 
specified in a basis of Killing vector fields. Since it is tensorial it is independent of basis and in terms 
of covariant differentiation along the world-line it may be written 

m ^A=-l-^t v (V A\/ V A). (45) 
o 47reo 

This is a system of third order differential equations for the worldline £ M (tt)- Regarded as an initial 
value problem it requires unfamiliar initial data (£^(0), £ M (0), £ M (0)) and solutions exist corresponding 
to self acceleration which must be regarded as unphysical. 

In the presence of an external Maxwell field jF ext one must confine the motion of the particle to a 
domain U that excludes the sources of T e *t- In this situation one expects that the equation of motion 
will acquire a contribution from the Lorentz force qoiv^ext' 

2 q 2 ~ 

m c 2 A = q i v Fe,t + r H V A (46) 
o 47reo 

so that in® is identified with the rest mass of the point particle. 

Although solutions to this system that self-accelerate can be eliminated by demanding contrived 
data at different points along the world-line there remain solutions that pre-accelerate in situa- 
tions where the external field is piecewise defined in spacetime Although the duration of pre- 
acceleration is probably classically unobservable for most electrons this general feature suggests that 
not all of the assumptions above are acceptable. 

One resolution of these difficulties [Tj is to assume that the right hand side of (J4T)|) should be 
expanded as a series in q with leading term for A given by -^iy-Fect- Then the above assumptions 
are held to be accurate only to some order in q and (jjS)) is understood as: 



A = -^ivFe* - q 2 2 f° iv(V A V v A cxt ) + ■ ■ ■ (47) 
m <r 6m c z 47re 

where A ext = ^1^^. The system is now manifestly a second order system of evolution equations. 
Although this offers a workable scheme it is unclear what its limitations are in different types of 
external field. Furthermore in situations where one has to contemplate the radiation from a large 
number of accelerating high-energy particles in close proximity the neglect of higher order terms in 
the expansion may be suspect. 

Given the complexities and reservations associated with ()47|) compared with the model leading 
to ([TTJl the latter is adopted in this article as a description of a collection of potentially radiating 
particles in the presence of an external electromagnetic field. Thus, our fundamental set of equations 
is 

dT = , (48) 

= -SlN*V , (49) 
eo 

V V V= — 2 iyT, (50) 

V-V = -l (51) 
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where X ■ Y = g(X, Y) for all vector fields X, Y e YTM. 

To simplify the subsequent analysis the system (J48H5ij) is recast using the 2-form F = -^sT and 

2 

the scalar field p = eo ^° e a -A/* where qoAf is the proper charge density since Af is the proper number 
density. This yields the system of field equations 

dF = , (52) 

d*F = -p*V, (53) 

V V V= i v F, (54) 

V-V = -\ (55) 

for the triple (V,p,F). Equation (J53j) leads immediately to the integrability condition (conservation 
of electric charge) 

d*(pV) = 0. (56) 

The scalar field p will be referred to as the reduced proper charge density. Since M is a number 
density, J\f > and so p > for q > and go < 0. 



2 Symmetric solutions to the coupled field system 

Although the field system (j52H55j) appears considerably less complicated than a large number of 
ordinary differential equations for a collection of accelerating charges, it is non-linear in (V, p, F) 
and so obtaining exact solutions is difficult except in certain simple circumstances. One method of 
obtaining solutions in more general situations is to employ a perturbation scheme based on properties 
of a particular class of exact solutions. 

In this section a family of exact dynamical solutions with symmetries generated by -J^ and -J^ will 
be described. Such highly symmetric solutions can be interpreted as accelerating "walls of charge" 
moving in vacuo under the influence of their self-fields (their "space-charge" fields) and an externally 
applied field. The coordinate y° will be identified with time t in an inertial laboratory frame (the 
speed of light c = 1) and x = y 1 , y = y 2 , and z = y 3 . The solutions will be independent of y 1 and y 2 
i.e. 

Cj^F = £ e F = 0, 

dy 1 9y' 2 

£ a V = £ a V = 0. 

dy 1 By' 2 

Thus one may reduce the problem to a field theory on a 2-dimensional Lorentzian spacetime with 
global coordinates y°,y 3 . The field system (|52H55|) is solved exactly using a co-moving coordinate 
system (r, er) adapted to the charged continuum. However, expressing the solutions in terms of 
laboratory coordinates (t,z) requires the inverse of the mapping (r, ex) — ► (t,z), which is generally 
difficult to obtain in closed form. To progress, a running parameter e > is introduced into the 
mapping (r, a) — ► (t, z) and a perturbation scheme in e is established. This approach facilitates 
an order-by-order construction of the inverse of the mapping (r, a) — ► (t, z) and thereby leads to 
1-parameter families (V £ , p £ , F e ) of solutions in e. It will be shown that 

oo oo oo 

n= — 1 n=— 1 n=l 
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over some range of e where the coefficients F n , V n and p n are 2- forms, vector fields and scalar fields 
respectively. 

Exact solutions to the system of equations (|52H55j) are obtained using the ansatz 

F = £(t,z) dtAdz, 

V = 1 (d t + /i(t, z)d z ) 

y/l-fl 2 (t,Z) 

where /i is the magnitude of the Newtonian velocity field of the charged continuum measured by the 
inertial (laboratory) observer d t . Equations (J52*)l and are manifestly satisfied while and (|51|i 
lead to 

dS = p#V, (57) 
y v V = £#V (58) 

where # is the Hodge map associated with the volume 2-form #1 = dtAdz. Since iyj^a = #(aAV) 
for all forms a independent of dx and dy, the action of iy on (|5T|) yields 

V£ = (59) 

i.e. £ is constant along the integral curves of V. Let C a be a 1-parameter family of proper-time- 
parametrised integral curves of V where each value of a corresponds to an integral curve of V: 

r -> (t = i(r,a),z = z(r,a)) 

and 

dt Oz 
V = C CT *<9 T = — d t + —d z 
or or 

where C aif is the push-forward (tangent) map associated with C a . Equation (JoT?|) is 

(C a *d T )£ = d T (C* a £) 
= 

where C* is the pull-back map associated with C a . Hence C*£ depends only on a and can be 
expressed using the function £ : R — > R where 

C(cr) = C*£ 

% \ (60) 

= £ {t(T,a),z{T,a)) . 

Using Vd t = Vd z = the (t, z) components of (JSHJ) are 

d 2 i , , , .dz, . 
— (r,a) = aa) — (r,a), 

j£{r,a) = (61) 
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whose particular solution satisfying the initial conditions 

i(0,a) = 0, z{0, a) = a, 

dt , , dz , 
-(0,^ = 1, ^(0,a) = 0, 

at r = is 

t = t(r, cr) = sinh(C(a)r) , (62) 
z = z(t, a) = —r^ [cosh(Cfcr)r) - 1] + a (63) 

where the function ( : R — > R must be supplied as data (i.e. ( is an ingredient in the initial 
conditions). Note that the charged continuum is taken to be initially at rest in the laboratory, i.e. 

V \ T=0 = (CU<9 T ) | T=0 = d t . 

More general initial conditions could have been used, but the above scenario is sufficient. 

After ( has been specified, expressions for V and F in the (t, z) coordinate system are obtained 
by inverting (jo^j) and (JBlIj) to give (r, a) in terms of (t, z): 

t = f(t,z), 
a = cr(t, z) 

and using (J6T1|) it follows that 

£(t,z)=((a(t,z)). 

An initial condition on £ is expressed using (. Equation (jo^|) is used to show that the local 
hypersurfaces t = and r = are equal and so (JBlIj) yields 

o"(0, z) = z. 

Thus, flnOJ) leads to 

£(0,z)=((a(0,z)) 
= C(z) 

on the space-like hypersurface t = 0. 

Not all choices for ( are admissible because, as noted in the derivation of (j52H55J) . the reduced 
proper charge density p is positive: 



p = -i v (#d£) > (64) 



using dHZJ) and (jSHD- At t = r = 



p(0, z) = -^#(^£(0, z)dt + «9 Z £(0, z)dz) 

= d z £(0,z) (65) 

= CO*) 
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where j^dt = —dz and j^dz = —dt have been used and ('(z) = ^( z )- Hence, employing (JfHJ) 

C'to > o 

i.e. ( is a monotonically increasing function. 

The coordinate chart (r, a) has important physical significance because it labels the streamlines 
of the charged continuum (the integral curves of V). Physically, the coordinate transformation 
(r, a) — > (t, z) is expected to be valid over the entire (r, a) plane. This will now be shown. 

The determinant A of the Jacobian matrix of the coordinate transformation (r, a) — ► (t, z) is 



A = det 



at as 

dt dz 

dcr dcr 



c(*) s 



+ cosh(^(cr)r) 



(66) 



Since C'(°") > and cosh(£(cr)r) > 1 it follows that 

A > 1. 

Monotonicity of ( implies that ( must vanish somewhere. Let Co G R be the point where C( cr o) = 
and note 

lim A = 1 + kVo)r 2 

(T— »<Jo L 

using (|66)1. Therefore, A is finite for all (r, a) G M 2 and the coordinate transformation (r, a) — > (t, 2) 
is valid over the entire (r, cr) plane. 
The implicit equation 

1 



a = z 



i + [c(a)n 2 -i 



(67) 



for a = a(t, z) arises by combining (jo^j) with (JBUj) . 

The initial electric field £(0,z) = ((z) = ( ext + Cseif(^) where Cext models an externally applied 
field and Cseif is the self-field of the charged continuum. The field C ext is constant because it satisfies 
the vacuum Maxwell equation 

dCcxt = 

inside the charged continuum (consider equation (joTj) with p = 0). The self-field £seif is to be defined 
in terms of the initial reduced proper charge density p(0, z): 



Cself(^) = 7^ 

and satisfies the asymptotic conditions 



/Z POO 
p(0, s) ds — I p(0, s) ds 
-00 J z 



(6* 



limCseif(^) = ~Q, 
lim Cseif(^) = ~\q 

z— >— 00 Z 



where Q = p(0, z)dz is the total charge per unit area of the continuum. 
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Figure 1: The time history of the initial charge distribution ()69|) with zero applied external field. 



Illustrations of the above results are shown in figures in figures Q an d 121 the external field is 
set to zero so the charge distribution p evolves entirely under its own self-field. In figures El and 0] 
the external field is set to a finite value. Figures El and H] show the streamlines of V in spacetime for 
each case. The initial charge distribution is taken as the Gaussian 



In the first case (figures Q and |2J) the initial Gaussian "wall of charge" evolves into two mirror-image 
"walls of charge" that propagate with equal Newtonian speed in opposite directions and tend towards 
the speed of light as t — > oo. In the second case (figures El and 0} the initial Gaussian "wall of charge" 
is accelerated by a combination of the external field and its self-field and approximately maintains 
its shape. 



p(0,z) = exp(-^ 2 ) 



(69) 



and, using (|55|). the associated electric self- field is 



Cself(^) 




The external field C ext is chosen as: 
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Figure 2: The lines of constant r and o generated by the initial charge distribution ()69|) with zero 
applied external field. Lines of constant r, except r = 0, are dark grey and lines of constant a (the 
streamlines of the flow), except a = 0, are light grey. The black line r = coincides with the £-axis 
and the black line a = coincides with the t-axis. 




Figure 3: The time history of the initial charge distribution ()69|) with finite applied external field. 
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Figure 4: The lines of constant r and a generated by the initial charge distribution with finite 
applied external field. Lines of constant r, except r = 0, are dark grey and lines of constant a (the 
streamlines of the flow), except a = 0, are light grey. The black line r = coincides with the z-axis 
and the black curve lying in the region z > is o = 0. 



2.1 Laboratory frame description 

Although the solutions above are exact in terms of (r, a), solving (|67J) for a(t, z) in closed form is 
impossible if ( is arbitrary. One way of tackling this problem is to introduce a running parameter 
e > 0, replace a by a 1-parameter family of functions a £ and solve (}6Tj) for a £ order-by-order in e. 

To effect the coordinate transformation from (r, a) to (t, z) ( is replaced with the 1-parameter 
family of functions ( £ : 



where Cext = ~C-i an d Cseif = Co- One may interpret (ffHjl by saying that external field effects (the 
first term -<C-i) dominate the effects of space-charge (the second term (o). The coefficients of the e 
terms in ()70|) have been labelled by the e-order of the term in which they appear; this convention 
will be used consistently throughout this article. 

A solution to (jHTj) can be obtained order-by-order in e using the series ansatz for t > 



The electric field component S £ (t,z) and 4- velocity field V s can be written entirely in terms of 



C £ (a) = -C-i + Co(a). 



(70) 



oo 




(71) 



n=0 



cr £ , t and z: 



S £ (t,z) = ( £ (& £ (t,z)) 



V £ = v/l + [C (a £ (t,z))]H*d t + C (cr £ (t, z))td z . 
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Using ([70)l. (J57j) . (|7T|). the choice 4 C-i > and t > 0, by equating orders in e it follows that 



a 6 (t,z) =z-t + 



e_ _ 1 + 2t( {z - t) 2 /v , m 
C-i 2tC»! 



-e 2 + 0{e" 



£%t, z) = -C-i + (o(z-t) + C ° ( * t] e + 0(e 2 ) 
£ C-i 



-C-i + Co(^ - 1) 



t{d t + d z ) + 



l + 2t 2 (>(z-t) a , Co(^-*)* 



2tC- 



C- 



(72) 
(73) 
(74) 



— ^2). Although V s ■ V £ = — 1 the (i, 2) components of diverge in the limits t — > 00 



where - d2 

and t — > 0. Divergences as t — > 00 are expected because the dominant term in the electric field is 
constant, and so the charged continuum is undergoing constant acceleration to leading order in e. 
Although V s diverges as t — > 00, the electric current J e = p e V e is 



J e = p ?V e = #d£ £ 

= C (z - t) (d t + d z ) + e^l^zA {dt + dz) + o(e 2 ), 



(75) 



where (|57j) has been used. Since J e = p e V^ e , J e = O(s ) and V^ e = 0(e~ l ) it follows that the reduced 
proper charge density p £ = 0(e). The leading order contribution to p e is obtained by comparing (J75j) 
with (|74*j) and using J = piV_i where 



V 



n=0 



It follows that 5 



n=-l 

«0 - «) 



n=l 



C-i* 



+ 0( £ 2 ). 



(76) 



(77) 



Equation (J75j) shows that J e is bounded to 0(e 2 ) as t — >■ 00. The divergence as t — > of the 
coefficient of £ in equation (|74|) stems from the square-root term in (|67j). Evidently, the domain on 
which the approximate solution is valid (the half-plane t > 0) is a subset of the domain of the exact 
solution (the entire (t, z) plane). This is not too surprising since the charged continuum is at rest at 
t — and the leading term J in the series for J e is light-like: 



Jo ' ^0 



0. 



There is no observer frame in which Jo is instantaneously at rest. 
This simple example suggests that the series 



J2 £np ^ 



ye 



J2 ^ 



n=-l 



;i = -l 



E 

n=l 



£ n pn 



where F_i is an external field (a solution to the source- free Maxwell equations), should be inserted 
into the field system (j52H55|) which is then solved order- by-order in e. 



4 Similar solutions are obtained for £_j < but with z + 1 replacing z — t. 

5 It should be noted that the equation p £ = — J £ ■ V e cannot be used to calculate (|77|l without first calculating J2 
since pi = -(J 2 • V_i + Ji ■ V + J ■ V\). 
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This approach can lead to the series expansions for the highly symmetric solutions just discussed, 
but can also yield e expansions for solutions that depend on (x,y) as well as (t,z). A consequence 
of the above series expansions is that the electric 4-current has the form 

oo 

Js = p eye = Y^e n J n . 

n=0 

It will be shown below that the above expansions partially decouple the field system (J52H55J) yielding 
an infinite hierarchy of equations that are amenable to solution when supplemented with appropriate 
boundary conditions and initial data. 



3 The perturbation scheme 



Based on the above exact symmetric solution it is assumed that there exists a class of solutions 
to (j52H55j) that represents configurations of charged particles in ultra-relativistic collective motion. 
It is shown in Appendix IA.1I that the notion of an ultra-relativistic 4-velocity vector field can be 
made precise and should be distinguished from the magnitude of a relative Newtonian speed that 
can be defined for any two time-like 4-vector fields. An ultra-relativistic velocity field is a pointwise 
limiting concept that depends on the existence of the forward light-cone structure at each event 
in spacetime. To leading order in the expansion defined below the velocity field of the charged 
continuum is light-like. The full series will be considered as an asymptotic expansion for a solution 
to (|52H55j) and physically represents an ultra-relativistic configuration controlled by the parameter 
e. Such configurations cannot of course be exhaustive. They are chosen to be representative of the 
class relevant to charged beams in high-energy accelerators. 
Introduce perturbation series for (V s , p £ , F £ ) in e of the form 



E 



e Pr, 



E 



(7f 



where 

V n G TTM , Pn G TA°M , F n G TA 2 M 
Theorem 1. Using \5 6 35b\) the coefficients of the expansions in fl7&| ) satisfy 

rfF n _ 1 = for n G {0, 1,2,...} 

{0 for n = 

_ 
r=l 



y~Vp r y n - r -i for n G {1,2,...} 

n 

fornG {0,1,2,...} 



r=0 
n 



r=0 



■1 for n = 2 



]^ Vr ~ 1 ' Vnrl= 10 fornG {0,1,3,4,5,...} 

n+l 

E d * (prK-r) = for n G {0, 1,2,...} 



(79) 
(80) 



(82) 
(83) 
(84) 



r=l 
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The partially decoupled equations \8O\84\) are amenable to an ordered analysis and can be arranged 
into the hierarchy: 



Solve flU and (JEU) with n = N for F n _ v . , 

Solve (El and (El with n = iV for K-i- 



For JV G {0,1,2,...} 
Step 3N + 1 
Step 3N + 2 
Step 3iV + 3: Solve (JEU) with n = N for p n+1 . 

Proof. Equations (j80H84j) are obtained by inserting the series (|75jl into (j52H55j) and the integrability 
condition (|56|) . The hierarchical structure follows by inspection. □ 

To illustrate the partially decoupled structure in (JHo"|) consider the first 8 steps in detail: 

Step 1: Adopt an external electromagnetic field F_i i.e. a solution of the source-free equations 

= and d*F^ = 0. (86) 

Step 2: Solve 

VviiV-i = iv-iF-i subject to V-i ■ V-i = (87) 
for V-i ^ 0, where F_i is data obtained in the previous step. 

Although (J87|) is non-linear in the unknown field V_i it is straightforward to analyse; it can 
be written as a quasi-linear second-order ordinary differential equation for the integral curves 
of V-i and is the only non-linear differential equation in this scheme. 

It may seem that the pair of equations above offer 5 scalar equations for 4 unknown scalars. 
However, the contraction of the first equation with V-\ identically vanishes, so there are only 
4 independent equations. This can be seen 6 by writing ()87|) with respect to an arbitrary 
basis and noting that one of the 5 components of (|HTj) is a consequence of the other 4. 

Step 3: The leading order reduced proper charge density pi is a solution to 



d* ipiV^J = 0. (88) 

Step 4: The 2-form F is a solution to Maxwell's equations with the current p±V^i as a source: 

dF = and d * F = - * Pl V- X . (89) 

The conservation equation (J88)) ensures that (|89|) has a solution. 
Step 5: The zero order velocity field Vq is obtained from the linear equations 

Vy_ i y + Vy V"_i = iy^Fo + i Vo F-i subject to V-i ■ V = 0. (90) 

Step 6: The second order coefficient p2 in the e expansion for the reduced proper charge density p e 
satisfies 

d * ( p 2 V^ +d* (p 1 V ) = 0. (91) 



6 The split of (|87|l with respect to the laboratory frame is given in Appendix IA.2I The equation obtained by 
contracting l|87|l with the laboratory observer field dt is redundant because it follows as a consequence of (|119|l and 
1(120)1 . 
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Step 7: Equation (|9*T|) ensures that 



dF l = and d * F a = - * p 2 ^i - *Pi Vb (92) 

may be solved for the coefficient F x (the first order term in the £ expansion of F £ ) given 
appropriate initial and boundary conditions. 

Step 8: The linear partial differential equation for V\ is 

Vv^Vi + V Vo V + V Vl V-! = i v _ 1 F l + i Vo F + i Vl F_ x subject to 2V^ ■ V x + V ■ V = -1. 

(93) 

l/ e must be calculated to at least first order in e to obtain a time-like vector field. To see this 
let X s = e^V^i + V and note X s ■ X s = V ■ V using the metric-product conditions in (J57)l and 
(I9UI) . Using lemma in Appendix IA. II Vn ■ Vq > and so X e is light-like or space-like. Now introduce 
Y £ = X £ + eV x = e~ V_i + V + and note Y £ ■ Y £ = V ■ V + 2V X ■ V- X + 0(e) = -1 + 0(e) using 
()93|) . Therefore, there exists some value of e > for which Y £ is time-like. 

It has already been stated that the pair of equations (|87|) in step 121 may be considered as 4 
non-linear equations in 4 unknowns. Furthermore: 

Lemma 2. TTie equations in each step 3N + 2, /or N > 1, are equivalent to a set of 4 independent 
inhomogeneous linear equations for the 4 unknown components o/V/v-i- 

Proof. Step 3iV + 2 involves solving (jg2*j) and (J53)l for Vn-i- Assume that the equations in steps 1 to 
3N + 1 have been solved and define L n as 

n 

Ln = J2 (Wn-^K-l - IVr-iFn-r-l) • (94) 
r=0 

Clearly (J52*|) is equivalent to L n = 0. The vector fields VLi, . . . , Vjv-2 have been obtained in steps 1 
to 3A" + 1 and so L n = for n — 0, . . . , N — 1. It will now be shown that given any Vn_i satisfying 
(J83J) with n = N the scalar field iy_ Y L N identically vanishes. 
From lj55j) for n = 0. . . . . N 

(n \ n n 

Y V r -l ■ V n -r-i = ^(Vy s K-l) • K-r-1 + Y K-l ' Vy s K-r-l 
r=0 / r=0 r=0 

n 

= 5^(Vv.v;-i) • K-r-i + 5^ • Vv.y g _i 

r=0 g=n 

where the second sum has been re-labelled by q = n — r. The two sums are identical so 

n 

Y(Vv s Vn-r-l) ■ Vr-l = 0. (95) 
r=0 
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Using fl§U) and L n = for all < n < N 



N 

n=0 
N N-n 



y^jv„-l ^ (Vv w _ n _ r _ 1 V r r _ 1 - iVr-± F N-n-r- 
n=0 r=0 
N N-n 

Z2 ^2 (Y n - 1 ' W N _ n _ r _iK-l - iVn-JVr-iFN-n-r-l) 
n=0 r=0 
N N 

Yl Yl t^- 1 ' V Wl^-n-l - iVn-Jv.-n-^N-a-l) (96) 



n=0 s=n 



where, in the final step, the innermost sum has been re-labelled by s = n + r. The double summation 
in (jnHJ) is over the set 

{(n,s) \n,s G {0,l,...,7V},n < s] 



and so 

N N N s 



ns 

n=0 s=n s=0 n=0 



YY Sns S '' 

n=0 s=r. 

for any summand S ns . Therefore, using (jnfij) 

N / s 

iV- X L N = I X/ ' ^jv-s-iVs-n-l - ^ ^V„-i^V a _ w -i-Pjy-a-l 

s=0 \n=0 n=0 

= (97) 

where the first term in the summand vanishes because of (J$H|) and the second term vanishes because 
ivn-Jvs-n-! = —ivs-n-i^Vn-i an d the sum over s is from to n. 

Note that the lemma applies to step 5 and beyond and so, for the purposes of the lemma, the fields 
F-i, V-i, pi, Fo obtained in steps 1 to 4 are prescribed data. Let {Xi, X 2 , X 3 , X 4 }, with X{ G TTA4 
for i — 1,2, 3, 4, be any vector frame such that X\ = V-i and define {e 1 , e 2 , e 3 , e 4 }, where e l G TA 1 ^, 
to be the coframe dual to {Xi, X 2 , X 3 , X 4 }: 

e*^) = 5) 

where <5] is the Kronecker delta. Equation (J%2~|) can be written as 

C(y N -x) = &v-i (98) 
where £ is the R-linear differential operator on TA 1 ^ 

where £at_i G TA 1 ^ depends on {V-i, Vq . . . , Vjv-2, F , . . . , Fjv-i} but is independent of Vjv-i- 
Hence, using (|9*7|) 

z Xl (£(^v_ 1 ) - ^-i) = WA^-i) - &r_i) = (99) 
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Likewise, equation (|83[) can be expressed as 

1 Xi Vn-1 = W-iVn-! = CV-i (100) 

where the scalar field C^-i € TA°A^ depends only on {Vo, . . . , Vat_ 2 }. 
Using the coframe {e 1 , e 2 , e 3 , e 4 } 

4 

Vjv-i = Ajv-i,je l 
for component fields Ajv-i,». From (jlOOj) Ajv-i,i = Cn-i and (|HH|) becomes 

4 

c(^2\ N -ij/) =£ N -i-C(C N . 1 e 1 ). (101) 

i=2 

From (j^nj) ixi annihilates ()101f) identically and 

^ A *-w e< )) = ix ' ^ N ~ x ~ c ( Cn -^) for * = 2 ' 3 ' 4 - 

Thus, using the inhomogeneous linear algebraic constraint (|83|). 3 inhomogeneous linear differential 
equations have been obtained for the 3 unknown fields {Aat_i )2 , Ajv-i,3, Ajv-1,4} and Ajv-i,i is obtained 
using (JHHJ), determining Vjv-i- □ 

In conclusion, the only non-linear equations in the hierarchy developed in theorem ^ are those 
for V-\. All of the remaining equations are linear and therefore, in principle, simpler to solve. 



4 Example : Charged beam propagating in free space 

As a simple application of the hierarchy (|80H84|) consider a high-energy charged beam propagating 
in free space. Unlike the exact "wall of charge" solutions derived earlier, the charge distributions 
considered here have finite extent transverse to and along the direction of propagation. The first 
eight steps of the hierarchy (|8UH84|) are: 

Step 1: The leading order Maxwell field F_i (the external field) is chosen to vanish so (|5B|) is trivially 
satisfied. 

Step 2: Using (|57jl the leading order 4-velocity field V_i satisfies 

Vv^V-i = 0, V-! • V-i = (102) 

since the external field F_i = 0. Therefore, consider solutions to (|8UH84j) adapted to a null 
geodesic coordinate system (u, v, x, y) where u = z — t, v = z + t and 

g = —dt (gi dt + dx <g> dx + dy <g> dy + dz <g> 
= - (dw <g> c?t> + (if <g> du) + dx ® dx + dy ® dy (103) 
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with the spacetime volume 4-form 

★1 = dt A dz A 

where 

= dx A 

The pair (dt, d z ) and their duals (dt, dz) are related to (d v , d u ) and ((if, dw) as follows: 

d v = - (d z + d t ), d u = ^ (d z - d t ) , 
dv = dz + dt, du = dz — dt. 

Metric identities used in the subsequent steps are 

*(du A a) = (-Ifdu A #j_a, *(cfo A a) = -(-Ifdv A # ± a (104) 

where #j_(a A X) = and a is any p-form such that z^a = ^ct = 0. The vector field 

V-i is chosen as 

V.! = 7-1 + d z ) 

= 2 7 _i3 t , (105) 

where 7_i 7^ is assumed and, using V<9„ = V<9 U = 0, the only non-zero component of ()102j) 
is 

9,7-i = 0. (106) 

Hence, 7 _i is independent of v. 
Step 3: Using (|105|) . equation (fSHJ) is 

d v (piy-i) = 

and so, using fjl(JCj|) 

^pi = 0. 

It follows that pi is independent of v. 
Step 4: Using ()104|) . particular solutions to (jBHJ) may be written 

F = d% A du (107) 

where the scalar potential $ is chosen to satisfy 

d v % = 0, d ± # ± d ± $ Q = Pi7-i#±l (108) 

and dj_ is the exterior derivative in the (x, y) plane defined by 

d±f = d x f dx + d y f dy 

for any 0-form / E TA°M. 
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Step 5: Using (|5Uj). the equations for V are 

Vv_ x Vb + Vv VLi = iv^Fo, ^-i-Vo = 
and are solved by the particular solution 

^o = 

since i v ^F = using (|T7T7jl . IjTTTEjl and (TTU51) . 

Step 6: Since Vo = it follows that has the form 

d* (p2^-i) = 
and, using (|1(J5|) and (jl(J6|) . is solved by any p 2 satisfying 

d v p 2 = 

i.e. p 2 is independent of v. 
Step 7: Since Vo = the Maxwell equations are 

dFi = 0, d * F 1 = -p 2 * Vli 

and, employing the same method used in step HJ are solved by the particular solution 

Fi = #iA(f« (109) 

where $i is chosen to satisfy 

d„$ 1 = 0, d ± # ± d ± $x = p 2 7-i#±l- (HO) 
Step 8: Using the equations for Vi are 

Vv_ i yi + V Vl VL 1 = iv_ 1 Fi, Vli-Vl^-^ (111) 

and since iy_ 1 Fi = (which follows using (jlU9j) . (|110|) and ()105|l ) a particular solution to 
(fTTTjl is 

Vi = 9 tt , 9 u 7_i = 

27-1 

using ([106)1 . ()105|) and V<9„ = V<9 U = 0. Hence, the function j_i is independent of u and 
therefore, using (|106jl . 7_i depends only on x and 

In summary these solutions, describing an ultra-relativistic charged distribution propagating along 
the z-axis with its electromagnetic self-fields, are 

V e = -2 7 _A - e-^—d u + 0(e 2 ) 
e 27_i 

^ + i^) s ' + (^-4^)* + °( £2 ) 

F £ = d% Adu + ed^> 1 A rfu + 0(e 2 ) 
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where <3> and $1 satisfy 



d±#±d±$o = 7-ipi#±l, d±#±d±$ 1 = 7_ip 2 #±l 



(112) 



and 



$ = $oK x, y), $1 = x, y), 

pi = pi(u,x,y), p 2 = p 2 (u,x,y), 
7-i = 7-i(x,y) 

with u = z — t. The scalar fields po,Pi,7-i are determined by their values on the space-like hyper- 
surface t = given as data. The potentials $o and $i are then solved in terms of 7_ipo and 7_ipi 
using the 2-dimensional Poisson equations (|112|) in the (x, y) plane. 

The 3- velocity of the beam is along the direction z in the laboratory frame with Newtonian speed 



h-x-e 1 



4 7 _i 



0(e 2 



i - 



27^i 



0(e 3 ). 



For example, consider a Gaussian bunch with transverse radius Rq travelling at constant Newtonian 
speed 1 — &x to order e 2 : 



l-i(x,y) = b , pi(z,x,y) = a exp 



x 2 + y 2 \ „ 



where ao, -Ro and 6o are constants and H : M — > K. is a smooth bump function vanishing outside the 
interval (— Z\,z\) and H(z) = 1 for z G (—2:2,22) and 21 > ^2 > 0. Then the laboratory reduced 
charge density 7-1P1 for some range of t is 



7-iPi = 7-i(ac,2/)/5i(z - t,x,y) = a b exp 



x 2 + y 2 
R 2 



E(z-t). 



Working in the cylindrical polar coordinates (t,R,<f),z) where x = Rcos<ft and y = Rsm<ft, a cylin- 
drically symmetric solution to ()112|) well-behaved at R = is 

* o= {r oo6 °^[ i " exp ("S)w E<2_i) 



and the corresponding electromagnetic 2-form F is 



R 2 

Fn = Of) bn — — 

2R 



1 - ex P ( -~B2 

K 0/ J 



E(z — t) dR A (-dt + dz). 



The laboratory electric field is radial, the magnetic field is azimuthal and their magnitudes are equal 
and vanish outside of the support of E. 



5 Conclusion 

When analysing the dynamics of charged particle beams it is often fruitful to adopt a description 
based on classical fields rather than classical point particles. Pathologies (such as pre-acceleration) 
associated with radiating point particles are avoided by relying on field-theoretical notions. 
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A novel analysis of charged beam dynamics has been presented and a model of a freely propagating 
charged bunch discussed. The approach relies on an asymptotic series representation of solutions 
to self-consistent spacetime covariant field equations for (V, p, F) describing a charged continuum. 
The asymptotic series for V is based on the light-like vector field V-i leading to an ultra-relativistic 
approximation. The hierarchy of equations obtained are more amenable to analysis than the original 
non-linear field system and particular solutions have been presented. 

There are numerous avenues for the development of this work involving ultra-relativistic charged 
beams in the vicinity of beam pipes, RF cavities, spoilers, etc. leading to dynamical effects that are 
often described in terms of "wake-fields" [HI E| • This work will lead to a clearer understanding of 
radiation-reaction exhibited by continuum models of charged particle beams. 
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A Appendix 

A.l Definition of ultra-relativistic and light-like limited vectors 

The spacetime description of a classical point particle with mass m > invokes a parametrised 
curve (world-line) with a future-pointing time-like tangent vector. For affinely parametrised curves 
the tangent vector V (4- velocity) is normalised: V ■ V = — 1. A spacetime reference frame may be 
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associated with a time-like vector field U with U ■ U = — 1 and an "observer" in such a frame is 
modelled by an integral curve of U. The 4-momentum of the particle is defined to be p = mo V. 
At an event where an observer curve intersects the world-line of the particle one has the orthogonal 
decomposition: 

p = £U + P 

where U ■ P = 0. Relative to U, the scalar £ is the energy of the particle (in units with c = 1) and 
P is its 3-momentum. The Newtonian velocity v of the particle relative to U is 

P V 

v = l = -vU- u (113 > 

Since U ■ U = V ■ V = —1 implies £ 2 = P • P + m 2 , it is traditional to say that the particle is 
"relativistic" relative to U if £ 2 ~ P • P or equivalently v ■ v ~ 1. Clearly, such a notion depends 
on both V and U and the nearness of to unity. Furthermore, one may always find a frame in 
which P = and £ = itlq. However, given a particular frame U one may contemplate a 4-velocity V 
that is relativistic relative to it. To generalise this notion for any frame one may consider a family 
of 4- velocities V £ and invoke properties of V £ as e tends to some limit. 

Let e € (0, £ max ) = { e ^ IK. 1 < e < e max } be the running parameter introduced in the main body 
of the article. A 1-parameter time-like vector field (which is not necessarily normalised) W £ is 
ultra-relativistic if for any time-like vector field Z G TTM 

f\W e ■ Z\\ 

lim MTT , I, = oo (114) 
e^o \ \\W £ \\ ) y ' 



where ||X|| = y/—X ■ X for any time-like (or light-like) vector field X e TTM.. Note that in this 
definition Z is not a 1-parameter vector field. 

Lemma 3. Our definition of an ultra-relativistic vector field on spacetime is independent of 

the choice of the time-like vector field Z. 

Proof. Clearly ()114|) is independent of any rescaling of Z. Therefore let Z and Z be two g-normalised 
future pointing time-like vector fields (i.e. Z ■ Z = Z ■ Z = — 1) and define the vector field X\ where 

Xi (Z + (Z-Z)Z) 



{{Z ■ Z) 2 - 1)1/2 
Hence 

Z = -(Z ■ Z)Z + ({Z ■ Z) 2 -l) 1 ' 2 X l (115) 

and Xi ■ Z — 0, X\ ■ X\ — 1. Choose space-like unit vector fields X 2 , X 3 so that (Z, X\, X 2 , X 3 ) is a 
^f-orthonormal frame. It follows that 

W £ ■ W £ = -{W £ ■ Z) 2 + (W £ ■ Xx) 2 + (W £ ■ X 2 ) 2 + (W £ ■ X 3 ) 2 

which implies 

(W £ ■ X x ) 2 _ (W £ ■ Z) 2 ( (W £ ■ X 2 ) 2 (W £ ■ X 3 ) 2 



\W £ \\ 2 \\W £ \\ 2 \\W £ \\ 2 \\W £ \\ 2 
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and hence 



l"l<ff^. (11 6) 



\w £ \\ ~ \\w e 



The metric contraction of pi5j) with w^iu is 

Z-W £ , A „z-w £ 



\\w £ \\ ' iiir-n ' ' ■ mr- 

and so 



'Z ■ Z« = -*J£ + ((Z ■ Z)» - l)**^ 



which, using the triangle inequality for scalars, implies 



iz . z|^ < Jfjn + (( z . z)» - iH*'T' 



Rearranging the above equation gives 



\\W £ \\ ~ ||iy £ || ||iy e || 

and so, using ()116j) 

I* ' ^ > |Z • Z|^iP - ((Z • Z)« - D" 2 ^? = ™ f|Z . z\ - ( { z ■ zf - iy» 



\\W £ \\ - 1 1 \\W e \\ vv 7 ' II „ 

Since \Z ■ Z\ — {{Z ■ Z) 2 — l) 1 / 2 > 0, and Z and Z are independent of s, it follows that 

,. \Z-W £ \ |Z.W e | 

hm — > oo ==>- hm — ► oo. 

e^O \\W £ \\ £^0 \\W £ \\ 



□ 



The connection between our definition of an ultra-relativistic vector field and v-v ~ 1 is exhibited 
by introducing the 1-parameter family of Newtonian velocities v e of W £ with respect to the observer 
field U: 



W £ 

v £ = - U 

W £ ■ U 



where U ■ U = — 1. It follows that 



£ e _ / W £ \ f W £ \ W £ ■ W £ 2(W £ ■ U) _ W £ ■ W £ 

V ' V ~\W r U + )'\W £ ~U + ) ~ {W £ ■ U) 2 + W £ -U + ' ~ {W £ -U) 2+ ' 

Thus lrnw, v e ■ v £ = 1 <^ lmw, W £ ■ W £ /(W £ ■ U) 2 = - lim £ ^ (\\W £ \\ /\W £ ■ U\f = and 
hence lim £ ^o v £ ■ v e = 1 -<=>■ lim e ^o \W £ ■ U\/\\W £ \\ — ► oo (i.e. W £ is ultra-relativistic). 

Observe that in definition ()114j) lim £ ^o W £ may or may not exist. Given a 1-parameter vector 
field, W e , W £ is said to be light-like limited if >V° = lim e ^ W e exists and is light-like (i.e. W° ^ 
and ||W°|| = 0). The relationship between ultra-relativistic vector fields and light-like limited vector 
fields is given by the following lemma: 
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Lemma 4. IfW £ is light-like limited then W e is ultra-relativistic. IfW 6 is ultra-relativistic and the 
limit lim e ^ W £ exists and is nowhere then W e is light-like limited. 



Proof. If W £ is light-like limited then \im £ ^ \W £ ■ Z\ = | \im £ _> (}V e ) ■ Z\ = |W° • Z\ > since 
W° is light-like and Z is time-like. Furthermore, lim £ ^o(VV e • W £ ) = — ||W°|| 2 = and so W e is 
ultra-relativistic. 

If W £ is ultra-relativistic and the limit W° = lim £ ^o VV £ exists and is nowhere then (since W° ■ Z 
is finite) ||W°|| must vanish to ensure (jll4|) is obeyed. □ 

The 1-parameter family of velocity fields V s introduced in the main body of this article (see section 
EJ) is ultra-relativistic. Moreover, the vector fields V-i, V and V\ then have certain properties: 

Lemma 5. 



• The 1-parameter family of vector fields V e is ultra-relativistic if and only ifV-% does not vanish 
anywhere. 

• If V s is ultra-relativistic then Vq is nowhere time-like and V\ does not vanish anywhere. 
Proof. 

• If V-i 7^ then for any time-like vector Z, V-\ ■ Z ^ so 

\\m{\V e ■ Z\J \\V £ \\) = lim^ 1 !^! • Z\) = oo 

e-+0 e^O 

since \\V £ \\ = 1 and so V £ is ultra-relativistic. If there is some igM such that V-\(x) = 
then 

\im(\V £ (x) ■ Z{x)\l \V £ {x)\) = Um(e°|Vo(a:) • Z{x)\) ^ oo 

e— >0 e— >0 

so V s is not ultra-relativistic. 

• Assume Vo is time-like. Since V s is ultra-relativistic it follows V-i ^ and, furthermore, since 
Vq is time-like V-i ■ V ^ 0; however, this conclusion contradicts (J§Uj) which states V and V-i 
are orthogonal. Hence, Vq is not time-like (i.e. Vq is space- like, light-like or zero). 

Now since V-\ is light-like and Vq is not time-like then in order to satisfy 2V-i ■ V\ + Vq ■ Vq = —1 
from V\ does not vanish anywhere. 

□ 



A. 2 The perturbation hierarchy in Gibbs' 3-vector notation 

To facilitate a comparison with existing fluid descriptions of beam dynamics [TTH ITT] that employ 
Gibbs' 3-vector notation the equations (j52H56j) in section |21 are here transcribed into that language. 
By splitting with respect to an inertial observer field -dt where 

g = —c 2 dt <S> dt + dx ® dx + dy ® dy + dz ® dz 
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one obtains 

m c 2 

V • e = 7p, 

Qo 

__ , 1 1 9c 

V x 6 = — pp + — — , 

go c z or 

„ Ob 

Vx e + ^7 =0 ' 

at 

V-6 = 0, (117) 
^dp + / _p_ ^ \ = / + x ^ 

2 . P P 

-1 + —22 = _1 ' 

m o^(7P) + V • Op) = 

7<9t7 + — ■ V7 = -f^-p • e (118) 
mo rriQC z 

where a ■ b is the standard Euclidean dot product of two 3- vector fields a and b. Note that (jll8)) can 
be obtained by taking the scalar product of p with the differential equation for p and the algebraic 
constraint on 7 and p in ()117|) . Therefore, like the final equation in (jll7j) (charge conservation), 
equation ()118|) is redundant. However, the two differential equations are treated very differently in 
the subsequent asymptotic analysis. Equation (jll8|) is automatically maintained order-by-order by 
the algebraic constraint on 7 and p but charge conservation must be imposed order-by-order. Hence, 
(j!18|) is not included in the following. 

The scalar field p is the reduced proper charge density introduced in the main body of the article 
and the MKS fields (p, e, b) and 7 are related to (V, F) by 

7„ P 

v= 1 -d t + 



and 



c m c 

F = —dt AE+ — *(dt AM) 
itlqc m V / 

where (P, E, B) are the vector fields 



P=(p-i)d x + (p- j) d y + (p- k) d z , 
E = (e • t) d x + (e • j) d y + (e- k) d z , 
B = (6 • t) d x + (b ■ j) d y + (b- k) d z 

and (i,j,k) is an orthonormal inertial 3-vector frame. 

The expansions for the 1-parameter families e e , b e , 7 s , p e and p £ corresponding to (fTSj) are 

e e = e _1 e_ 1 + e + ee 1 + . . . , b e = e^b^ + b + eb x + . . . , 

7 e = £ -1 7-i + 7o + ^7i + • • • , V s = £~ l P-i + Po + epi + ■ ■ ■ , 
p £ = epi + e 2 p 2 + ... 

Substituting the above expansions into (J117)) yields the following hierarchy corresponding to equations 
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Step 1: The source-free Maxwell equations are satisfied by the prescribed external electric field e_i 
and external magnetic field b i. 

V • e_i = 0, 

1 <9e_! 



V x 6_i 



c 2 ^ ' 



V x e_i + = 0, 



V • 6 i = 0. 

Step 2: The partial differential equation 



+ f — ' V ) P-i = *> f 7-ie-i + — P-i x b_x ) , (119) 
dt \m J V m 



and algebraic condition 



Pi Pi 



Tii + ^ = (120) 



are solved for 7_ x and p i. 
Step 3: The leading order equation describing charge conservation is 



d 

m ^(7-iPi) + V • (pip-i) = 



and is solved for p\. 



Step 4: Solve the Maxwell equations for e and b with 7_ipi and ^p 1 p_ 1 as sources: 

m c 2 

V • e = 7-iPi, 

Qo 

, 1 1 <9e 

V x b = — PiP-i + ^^r> 

g c 2 at 

Vxe„ + f = 0, 

V • 6 = 0. 



Step 5: Solve 



dp dp-! f po \ ( p i 

7-i-aT + ^°^7~ + V p_i + V po 

at at \m I \ m 



go (7-ieo + 7oe-i) + — (P-i x 6 + p x b_ x ) , 

m 



p i • Po 

- 7-i7o + = 



ITLqC 2 



for 7q and p . Use of the above algebraic equation to eliminate 7 leads to an inhomogeneous 
linear partial differential equation for p . 
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Step 6: The equation 

d d 
mo — ^opi) + m ^(7-ip 2 ) + V • (pip ) + V • (p 2 P-i) = 

enforcing charge conservation is solved for p 2 . 
Step 7: The electric and magnetic fields e x and bi are solutions to the following Maxwell equations: 

_ m c 2 . 

V • d = (7-1P2 + 7oPi) , 

V X &! = — (P2P-I + PlPo) + 



go W ^- J ' ' r-> ■ 

V • 61 = 0. 



V x ei + = 0, 



Step 8: Solve 



dpi dp dp-! ( p_i \ / p _\ . / Pi _ , 



27-i7i - 7o7o + 



= g (7-1^1 + 7oe + Ti^-i) + — (P-i x f»i + Po x &o + Pi x b-i) , 
2p-i • Pi + Po • Po _ 

2 9 



for 7x and p 1 . Eliminating 7! using the above algebraic equation leads to an inhomogeneous 
linear partial differential equation for pi. 
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